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In this work we investigate the role of the beam astigmatism in the Goos-Hanchen 
and Imbert-Fedorov shift. As a case study, we consider a Gaussian beam focused 
by an astigmatic lens and we calculate explicitly the corrections to the standard 
formulas for beam shifts due to the astigmatism induced by the lens. Our results 
show that astigmatism may enhance the angular part of the shift. 


I. INTRODUCTION 

Geometrical optics considers light fields as rays directed along the propagation direction 
of the held itself. Within this approximation, most of the phenomena that we witness daily 
regarding light can be easily explained. An example is given by the phenomenon of rehection 
from an interface, which for an optical ray happens in a specular way, following Snell’s law [1]. 
When the wave properties of optical helds are taken into account, however, deviations from 
specular rehection can be observed. This is the case of optical beams, whose hnite transverse 
sizes ahect its rehection and refraction across interfaces. The most common manifestation of 
these ehects is given by the so-called Goos-Hanchen [2-4] and Imbert-Fedorov [5-15] shifts. 
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the former occurring in the plane of incidence and the latter in the plane perpendicular to 
the plane of incidence. These phenomena have been extensively studied in the past for a 
vast category of beam conhgnrations [16-20] and interfaces [21-24], A comprehensive review 
on beam shift phenomena can be fonnd in Ref. [25]. Recently, the analogy between beam 
shifts and the qnantnm mechanical weak measnrements has been also pointed out [26-30], 
which resulted in the possibility of observing amplihed beam shifts [31, 32]. 

In this paper we theoretically investigate the effect of the beam astigmatism on both Goos- 
Hanchen and Imbert-Fedorov shifts. Onr calculations show that the astigmatism affects, 
at the leading order, both the spatial and angnlar shifts, with its main action being the 
introdnction of an enhancement factor in the angular part of the shifts. The results presented 
here address the simple case of an astigmatic Gaussian beam, that we take as paradigmatic 
example for stndying the effect of astigmatism on the Goos-Hanchen and Imbert-Fedorov 
shifts. 


II. MODEL AND METHODS 

We start our analysis by considering a monochromatic fundamental Gaussian beam of 
freqnency u = ck (being k = 2 ti/\ the wave number) characterized by a beam waist Wq, 
which impinges upon an astigmatic lens with two different focal lengths, namely fi in the 
longitudinal {x-) direction and ft in the transverse {y-) direction, as sketched in Fig. 1. 
After the lens, the beam will be focused at different lengths in the x and y direction. To 
calculate the beam waist in these directions, we make use of the well-known formula for 
paraxial beams [36]: 

—— = —-- ( 1 ) 

q{s) // ^ ^ 

where y = {/, t}, is the focal length, s is the distance between the Gaussian beam’s waist 
and the lens, and q{s) = s — izr is the complex beam parameter, being zr = kwl/2 the 
Rayleigh range of the beam [36] and g^(s + 2 ;) = qf^{s)+z. Obtaining g^(s) from the previous 
equation and setting its real part to zero gives the distance from the lens to the new waist 
as 

^ ^ _ g - /m 


( 2 ) 
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while the imaginary part of allows us to write the longitudinal and transverse beam waist 
as follows: 

/^ ^ fl 

Uoy' + 

According to the convention used in Refs. [25, 33], we describe the reflection process 
by means of three different reference frames, namely one attached to the incident beam 
Zi}), one attached to the reflected beam {{xr^Vri Zr}) and a third reference frame 
{x,y,z}, called the laboratory frame, attached to the reflecting interface. These reference 
frames are linked together by a rotation of an angle 9 around the y axis [37]. The complex 
electric field impinging on the dielectric interface can be written in terms of its Fourier 
components, in the incident reference frame, as follows: 

E^(r) = J2[ (4) 

A=1 

where d^K = dUdV, 9 is the angle of incidence and U,V and W are the components of the 
wave vector in the incident reference frame, defined by the relation 

k = k{Uici + Vji + Wzi) = + kyf + k^z. (5) 

In Eq. (4), Ax{U, V, 9) = a\{U, V, 9)A(U, V)ex(U, V, 9) is the vector spectral amplitude of 
the field, with ex{U,V,9) being the local reference frame attached to the incident field [35], 
ax{U, V,9) = f ■ ex{U, V, 9) are the projections of the initial beam polarization f = fpk + /^y 
(with l/pp + l/sp = 1 ) onto the local basis ei(k), 62 (k) and k attached to the wave vector 
k, and A{U, V) is the scalar spectral amplitude associated with the beam, that in our case 
assumes the form of an astigmatic Gaussian beam, i.e., 

A{U, V) = ( 6 ) 

where wqi and w^t are defined according to Eq. (3). 

If we assume the beam to be well collimated, the paraxial approximation holds, and the 
expressions for the local reference frames and the polarization functions can be written in 
the following simple form [35]: 

ei(f/, V,9) =Xi + V cot 9yi - Uzi, 


62 ( 11 , V, 9) = -V cot 9xi + y* - Vzi, 


(7a) 

(7b) 



for the local basis, and 
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a,{U,V,e) = fp + fsV cote, 

(8a) 

a,{U,V,e) = fs-fpV cote, 

(8b) 


for the polarization functions. Notice, moreover, that the expression for A{U, V) as given 
by Eq. (6) has been written in correspondence of the beam waist after the lens. The beam, 
therefore, has propagated a distance D = d — s — si (see Fig. 1) before reaching the reflection 
surface. This propagation factor is correctly taken into account by the 2 :-dependent part of 
the angular spectrum of the incident held. 

Upon rehection, the single plane wave components of the held described by Eq. (4) 
experience geometrical rehection, according to the Snell’s law [1]. The rehected electric held 
can be then represented by Eq. (4) with the substitution 

e,(k)e**^'^ ^ r,(k)e,(k)e**^•^ (9) 

where rx(k) are the usual Fresnel coefficients for p (A = 1) and s (A = 2) polarization [1], 
and k = k — 2z(z ■ k) = — + Uy^ + fUz^. The electric held in the rehected frame can 

now be dehned as follows: 

E^(r) = ^ / d^KA.iU, (10) 

where X = kxr, Y = kyr, Z' = k{zr + D) and the vector angular spectrum in the rehected 
frame is given by Ax{U,V) = A{U,V)rx{U,V)a\{—U,V,7r — 6)ex{—U,V,7r — 6). 

The beam centroid can be then calculated as the weighted average of the position vector 
R = Xxr + Yyr with respect to the total beam intensity in the rehected reference frame 
|E^(a;, p, z)\‘^. By employing the quantum notation for optical beams developed in Ref. [38], 
its explicit expression reads 

= (eV^) " 

where {^) = with ^ = X,Y. The spatial (A) and the angular (0) Goos-Hanchen and 

Imbert-Fedorov shifts are then given by the following relations: 

Ag// = (X)|2=o, ©GH = 

A/f = (U)|^=0; Q/F = 1 ^ • 


( 12 ) 

(13) 
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By substituting into Eq. (11) the expression of the electric held in the rehected frame 
given by Eq. (10) with a hrst order accuracy in [U, V) for the numerator and a second 
order accuracy for the denominator, we arrive, after a lengthy but simple calculation, to the 
following results: 


^GH — [(hEp0p + hEj^s) + Ej; {WpPp + WsPs)] , 


Qgh = 


1 


kz 


R 


\ ^ 

— iWpPp + WsPs), 

Wol / 


(14a) 

(14b) 


for the Goos-Hanchen shifts, and 

cote 


Aif = 
Qif = 


k 


sin,^ + 


G(K - 

RpRs 




kz 


R 


Wot 


RpRs 


cot 9, 


(15a) 

(15b) 


for the Imbert-Fedorov shifts. In the previous equations, px = Re{(91nrA/i96*}, (px = 
Im{(91nrA/(96*}, Rx = |rA| (with A G {p, s}) and 

D 2Re{qp{d)} 






(16) 


Im{g^(ci)} ’ 

where qp{d) can be calculated from Eq. (1). Moreover, IEa = fx^x/ifp^l~k fgR^^ + S) is the 


fractional power contained in each polarization corrected by the beam’s astigmatism factor 
S, whose explicit expression is given by 


5 = 






Wot f 


d^Rl 


2wli V de^ 


-2Ri] +coRd{Ri-R 


+ cot ORp^ +/2 
+ cotORs^-Rl 


w, 


0 , a^R. 


2w‘li 89“^ 


+ cot^ 9 


R: 


(17) 


Equations (14) and (15) are the hrst main result of this work. 


III. RESULTS AND DISCUSSION 

A closer inspection on Eqs. (14) and (15) reveals that the ehect of the astigmatism 
on the Goos-Hanchen and Imbert-Fedorov shifts is twofold. Firstly, the fractional power 
normalization factor Wa is changed by the quantity <5. This means that, in contrast to the 
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case of a fundamental Gaussian beam [33], the astigmatism introduces a correction on the 
energy stored in each polarization, being this correction dependent on both polarizations, 
as it appears clear from Eq. (17). However, since this correction is of the second order into 
{U,V), it is very small and can be therefore neglected. By doing so, therefore, the spatial 
shifts corresponds exactly to the ones of a fundamental Gaussian beam described in Ref. 
[33], 

A second, and more interesting effect of the astigmatism can be observed in the angular 
shifts, where an extra multiplicative factor appears. By neglecting S in the expression of 
W\, the angular shifts can be therefore rewritten in the following inspiring form: 

Qoh =(—) (18a) 


0 /F=( —joj?, (18b) 

where &gh if rsfei' fo the angular Goos-Hanchen and Imbert-Fedorov shifts for a funda¬ 
mental Gaussian beam [33], respectively, and their explicit expression, in our case, is given 


QgH = T- iWpPp + WsPs) , 

KZr 


0(0) - 
yjiF — 


'VW, {Rl - Rf) 


RpRs 


(19b) 


This is our second main result. The astigmatism introduces an enhancement factor 




{s - fpf + zj 


in the angular shifts that depends essentially on the degree of astigmatism of the beam in 
each direction, which in this case is here represented by the longitudinal or transversal focal 
length, respectively. 

It is worth noticing, moreover, that the main effect of astigmatism is to introduce a 
selective amplihcation of the angular GH and IF shifts, as it appears clear from Eqs. (18). 
According to those equations, in fact, while the amplihcation factor in front of the angular 


GH shift only contains the longitudinal focal length through the quantity wqi, the transverse 
focal length only appears (through the term wot) in the angular IF shift. 

In order to quantify this enhancement, in Fig. 2 G; is plotted against the normalized focal 
length ( = fi/zR. As can be seen, for values of the focal length (either in the longitudinal or 
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in the transversal direction, depending on which shift is one interested in) approaches zero, 
the enhancement factor diverges, thus giving a signihcative enhancement of the angular shift. 
Although a rigorous vectorial analysis of the focusing process for / —)■ 0 would reveal that 
the considered divergence is only an artifact of the paraxial equation, the essential feature 
of it, namely the amplihcation of the angular shifts due to a shorter focusing length of the 
lens, is still very well captured by our paraxial model. 


IV. CONCLUSIONS 

In Conclusion, we have investigated the effect of a beam’s astigmatism on the Goos- 
Hanchen and Imbert-Fedorov shifts that occur when the beam impinges onto a dielectric 
surface. We have shown that the main effect of astigmatism is the introduction of an 
enhancement factor on the angular shifts that is essentially depending on the degree of 
astigmatism of the beam in the longitudinal (transverse) direction. We have also shown that 
for small values of the focal length (compared with the beam’s Rayleigh range zr) of the 
lens we used to model the beam’s astigmatism, a drastic amplihcation of the angular shifts 
occur. We believe that this result could be used experimentally to have another available 
channel of amplihcation, that could contribute to easily realize gigantic beam shifts that 
could lead to practical applications of this phenomenon. 
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FIG. 1: Scheme of the beam reflection at the plane interface. Here fi and ft are the longitudinal 
and transverse focal lengths of the astigmatic lens, respectively. Note that for a cylindrical lens 
one has either ft = oo or fi = oo. In addition, s is the distance from the waist of the incident beam 
to the lens, and d — s is the distance from the lens to the reflecting surface. si is the distance from 
the lens to the longitudinal waist wqi. Similarly, we can define as st the distance from the lens to 
the transverse waist WQt (not shown in the figure). Finally, Zr is the distance from the reflecting 
surface to the detector position. 
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FIG. 2; Plot of the longitudinal enhancement factor 11; as a function of the normalized focal 
distance ( = fi/zR. As it can be seen, as /; —>■ 0, the enhancement factor grows, thus providing a 
huge amplification of the angular shift. The same argument is valid for 11; as well, provided that the 
substitution /; —)• ft is understood. The inset shows that when —)• oo (i.e., when /^ —)■ oo, with 
/i = the enhancement factor tends to unity. These plots have been obtained assuming 

S/ZR = 1. 










